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Целью работы является рассмотрение проблемы устойчивости 
установившихся стационарных траекторий упругих деформа- 
ционных смещений инструмента при продольном фрезерова- 
нии концевыми фрезами. В отличие от ранее рассмотренных 
проблем устойчивости, в которых анализируется случай мед- 
ленных движений, для которых в уравнениях в вариациях 
относительно стационарной траектории параметры системы 
можно рассматривать замороженными, анализируется случай 
большой скорости резания. В этом случае при анализе устой- 
чивости необходимо рассматривать линеаризованную систему 
в вариациях с периодически изменяющимися коэффициента- 
ми. По мере увеличения скорости вращения инструмента во 
многих случаях наблюдается параметрическое самовозбужде- 
ние колебаний, поэтому главное внимание уделяется изуче- 
нию условий параметрического самовозбуждения динамиче- 
ской системы фрезерования концевыми фрезами. Показано, 
что на условие параметрического самовозбуждения оказыва- 
ют влияние как технологические режимы резания, прежде 
всего частота вращения инструмента, так и геометрия инстру- 
мента, которая определяет матрицы угловых коэффициентов 
ориентации сил резания. Приводятся примеры областей 
устойчивости в зависимости от изменения параметров систе- 
мы. 


Ключевые слова: процесс фрезерования концевыми фрезами, 
стационарные траектории, периодически изменяющиеся па- 
раметры, устойчивость, параметрическое самовозбуждение. 
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Введение. При изучении динамики процесса резания, в частности, фрезерования, рассматривают упругие 
подсистемы со стороны режущего инструмента и обрабатываемой детали, которые взаимодействуют между собой 
через динамическую связь, формируемую процессом обработки [1-9]. В свою очередь, динамическая связь характери- 
зует модель сил резания, представленную в координатах состояния системы и технологических режимах [10, 11]. В 
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зависимости от модели сил, зависящей от конкретного процесса и условий обработки, анализируются проблемы 
устойчивости и многообразий, формируемых в окрестностях равновесия. Рассматриваются автоколебания [12, 13], 
инвариантные торы [14—17] и хаотические аттракторы [16-19]. Показано, что для раскрытия основных свойств упру- 
гие деформационные смещения не могут быть скалярными. Из всех процессов наиболее сложным является процесс 
фрезерования, который является нестационарным с периодически изменяющимися параметрами длины и толщины 
срезаемого слоя каждым режущим лезвием фрезы [20—41]. Поэтому в подвижной системе координат, движение кото- 
рой определяется траекториями движения исполнительных элементов, стационарным установившимся состоянием 
является не точка равновесия, как это имеет место при точении, а некоторая замкнутая траектория. В связи с этим 
уравнение динамики, в общем случае, имеет периодически изменяющиеся коэффициенты. Ранее был рассмотрен слу- 
чай «медленных» стационарных траекторий [41], для которых в рассматриваемый момент времени параметры дина- 
мической связи можно считать замороженными. Это справедливо для малых скоростей резания и достаточных пара- 
метров жесткости инструмента, что является традиционным случаем обработки. Однако по мере увеличения частоты 
вращения инструмента и (или) уменьшения изгибной жесткости концевой фрезы, частоты собственных колебаний 
инструмента и его вращения становятся соизмеримыми. В этом случае в динамической системе резания возможно 
параметрическое самовозбуждение. Свойства системы в этом случае принципиально отличаются от свойств, рассмот- 
ренных ранее. Именно этому случаю посвящена данная работа, которая фактически является продолжением статьи, 
опубликованной ранее [41]. 

Уравнение в вариациях относительно стационарной траектории. При увеличении частоты вращения 
шпинделя наступает момент, когда время периодического изменения коэффициентов становится соизмеримым с им- 
пульсной реакцией системы. На практике такие условия соответствуют скоростному фрезерованию. В этом случае 
исходное дифференциальное уравнение несколько упрощается, но его анализ становится более сложным. Во-первых, 
при высоких скоростях резания запаздыванием вариаций сил по отношению к изменениям колебательных смещений 
инструмента относительно заготовки можно пренебречь, во-вторых, для анализа устойчивости при построении урав- 
нения в вариациях относительно стационарной траектории необходимо учитывать, что параметры системы становятся 
функциями времени. В этом случае необходимо учитывать параметрическое самовозбуждение системы. Параметри- 
ческие явления наступают и в том случае, если фреза обладает малой изгибной жесткостью. Исходное уравнение ди- 
намики в этом случае, в отличие от рассмотренного в [41], имеет следующий вид (рис.1): 
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где Х,, Х› — упругие деформационные смещения инструмента в плоскости, нормальной к оси его вращения в си- 
стеме координат, привязанной к вершине режущего зуба; т, й, с — инерционные, диссипативные и упругие коэффи- 
циенты, равные между собой, так как концевая фреза обладает свойством симметрии геометрических параметров; 
Ь — ширина срезаемого слоя зубом фрезы; ру — давление стружки на переднюю поверхность режущего лезвия зуба 


фрезы; ау (2) —ЩЩ заданная функция кинематической составляющей вариации толщины срезаемого слоя; 


ть 2, хо . ии == угловые коэффициенты ориентации главной составляющей силы резания и силы, формируемой в 


области контакта задней грани инструмента с заготовкой; О , > — параметры скоростной зависимости сил; Е® 


значение силы, формируемой в области задней грани инструмента при скорости движения инструмента в сторону за- 
готовки без учета скоростей упругих деформационных смещений. 

На приведенной иллюстрации справа жирной линией показана функция изменения кинематической толщины 
срезаемого слоя а, (ф), то есть его изменения без учета упругих деформационных смещений инструмента. Там же 
дана пунктиром функция изменения толщины срезаемого слоя с учетом упругих деформационных смещений для 
установившегося состояния а(ф). 
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Рис. 1. Схема формирования сил и толщины срезаемого слоя в зависимости от угла поворота фрезы 


В системе (1) модуль главной составляющей силы определяется выражением: 
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Е® = Е® -ехр[-9>4Х/ / 4]. При построении уравнения (1) принято также, что частота вращения инструмента явля- 


ется постоянной 62 = соп5!. Поэтому угол поворота инструмента ф отличается от времени коэффициентом. Пусть за- 
> * * *ж2Г 
дана стационарная установившаяся траектория деформационных смещений инструмента Х’={Х!,Х.} . Она вычис- 
ляется по методике, изложенной в [41]. Тогда функции ах (1), Х “(0), АХ" / 4 являются заданными. 
Для рассмотрения устойчивости траекторий Х“ ={Х!,Х> } необходимо исследовать уравнение в вариациях 
относительно этих траекторий [42]. Представим деформационные смещения в виде Х, =Х, (+х (1) 
м < 
иХ, =Х.,()+х2(0. Здесь м({), х›(Р) малые вариации упругих деформационных смещений относительно траекто- 


рий Х | (1), Х а (Г). В свою очередь, функции Х 1 (1), Х > (Г) — есть установившиеся стационарные траектории. Следова- 


тельно, для скоростей упругих деформационных смещений справедливо 
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где У(@),У,() — заданные скорости деформационных смещений для стационарных установившихся движений; 
4 4х - 
а. их вариации. Таким образом, подлежащее анализу при изучении устойчивости линеаризованное уравне- 
ГА 1 
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где т= ; И = р — } 
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В 1(9) = Е Фа, ехр-а5М, 5 (2); И, (6) = ХЕ Фа, ехр[- 05, 5 ()]; №1 (9) = Хировоца; ()ехр/--04 У, (1); 
й 2 (1) = ХоРовоца$ (В) ехр[-о4 Г, 5 (1)] ; с 1 (1) = ХРов + пехр[-ои 7, ; (0)]} ; 1.2 (0) = ХоРов 1+ нехр[-04 У, ;(1)]} ; 
а; =а(-Х! 0; У ;@=У-У, @); И; @ =у(Р —У, (@) , У“ — постоянная скорость подачи. В уравнении (3) 


отброшены члены второго порядка малости. Проанализируем переменные коэффициенты в (3). Коэффициен- 
ты И, 1 (1), №>1(0), № (1), № (1) ‚11 (1) ‚ с, 2 (1) состоят из постоянной и периодически изменяющейся составляющих. Их 


период Тоодин и тот же. Период не меняется и в мультипликативных формах #1(1), #й,››(1), так 
как а; (1))0, ехр[-оц Г, ; (1)])0. Поэтому их можно разложить в ряд Фурье, в котором присутствуют постоянная со- 


ставляющая и только синус-преобразование. Ряд достаточно быстро сходится, поэтому все коэффициенты можно 
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представить в виде постоянных и периодически изменяющихся коэффициентов, имеющих лишь коэффициенты при 
синусах. Таким образом, базовое уравнение, подлежащее анализу, имеет структуру 
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Анализ условий параметрического самовозбуждения. Вначале рассмотрим несколько частных случаев. Во- 
первых, при больших задних углах режущих лезвий фрезы силами, действующими на заднюю грань инструмента, 


можно пренебречь, то есть Ё °) —=0. Тогда из (3) имеем # ;(#) = >(#) =0. Во-вторых, коэффициенты Х1,Х› принци- 


пиально зависят от переднего угла зуба фрезы. Рассмотрим два предельных случая. Первый случай 7, =1.2 =0. 


Второй случай 7, =0, > =1. В первом случае с учетом ЕЁ © =0 уравнение (4) преобразуется к виду 
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В системе (5) второе уравнение не имеет координаты д; ‚ и оно характеризует асимптотически устойчивую си- 


стему для х› . Поэтому х. (1) 5 = 0, а траектория х; (1) становится аттрактором. Время, в течение которого коор- 


при: 
дината х› =0 , зависит от коэффициента й. Начиная с этого времени система (5) вырождается в следующее скалярное 
уравнение 
а» 
а? 
Его свойства будут определяться хорошо изученным уравнением Матье-Хилла [43]. 
Во втором случае (4) преобразуется к виду 


4х 4х 
а (7) 





Ей | +{с+с [1+ Эзт(ОБ р =0. (6) 


0 с 0 
где т = ; (В) = ; (= р 
р . № | п) 1+52 зико |” ° с + 351(09] с 
Асимптотически устойчивой будет координата х,. Следуя изложенному выше синергетическому принципу 


анализа, получаем уравнение, свойство которого определяет устойчивость системы (4), также определяется скаляр- 
ным уравнением 
а?х . Ах 
т +{п-1 1+5 зп} 2+ =0. (8) 
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Система (8) может стать абсолютно неустойчивой при условии вв 1+5 $1(О1)] при 1 Е[0О, о). Поэтому 


нас будут интересовать, прежде всего, случаи, когда на некоторых промежутках времени указанное выше условие не 
выполняется. 


Приведенные выше случаи являются скорее исключением, чем правилом. В общем случае )! 1 Х> 1 и 


1 = 0, Х2 =0, а между ними существует известное условие нормировки (х: а (2)? =1. При реальных значениях в 


системе (4) за счет взаимного влияния двух колебательных контуров формируется сложный динамический процесс, 
который стягивается по своим свойствам к двум процессам, описываемым уравнениями (6) и (8). Система (4) прин- 
ципиально отличается от уравнения Матье-Хила, во-первых, тем, что в системе (4) имеет место взаимодействие двух 
колебательных контуров. Во-вторых, в этой системе периодически меняется матрица скоростных коэффициентов. 
Наконец симметричная часть матрицы скоростных коэффициентов может на отдельных временных промежутках ста- 
новиться отрицательно определенной, то есть ускоряющей. Тогда на этих временных отрезках система теряет устой- 
чивость. Причем потеря устойчивости, как правило, является колебательной. Тогда за счет чередования устойчивых и 
неустойчивых отрезков времени в системе формируются периодические движения. В общем случае для исследования 
системы (4) можно воспользоваться теорией Флоке для Ту -периодических процессов [43]. 
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Проанализируем условия параметрического самовозбуждения на примере системы, параметры которой при- 
ведены в табл. 1. Дополним эти данные составляющими матрицы упругости и скоростных коэффициентов, изменяю- 
щимися по синусоидальному закону, причем период их изменения одинаков. Построение диаграмм параметрического 
самовозбуждения выполнено на основе прямого интегрирования рассматриваемых дифференциальных уравнений по 
методу Рунге-Кутта 4-го порядка. Погрешность метода не превышала 0,5%, что оценивалось на основе сравнения из- 
вестных решений с полученными методом Рунге-Кутта. Построение диаграмм выполнено по сетке. 

Таблица 1 





т„[кг-с? [мм] | №[кг-с/ мм] | с[ке / мм] о [с/ мм] и Е® ке] 


3.0.103 5,0 3000 65 0,5 2,0 





























Выполним анализ по следующей схеме: 1) проанализируем область устойчивости для случая 
и; ; = =0,1,5=1,2; 2) для параметров из области устойчивости рассмотрим диаграмму устойчивости для уравнения 
Матье-Хилла (6); 3) проанализируем диаграмму параметрического самовозбуждения для системы (8); 4) для системы 
(4) рассмотрим два случая (х, =0,6, 5 =0,8 и Хх, =0,8, Хх. =0,6); 5) выясним влияние на параметрическое самовоз- 
буждение динамической связи, формируемой контактом задней грани режущего лезвия фрезы. 

На рис. 2, а приведена область устойчивости системы в плоскости двух параметров: ру и © системы (3). Ши- 
рина срезаемого слоя Б = 20мм . Слева от фигуративной линии система является устойчивой. Область вычислена для 
случая, когда и;, =3=0,1,5 =1,2. В дальнейшем рассмотрим две точки («1» и «2» на рис. 2), которые находятся в 


области параметров, при которых система устойчива. Однако точка «1» находится в непосредственной близости от 
фигуративной линии. На этой же иллюстрации (верхняя иллюстрация на рис. 2, 6 приведены диаграммы параметри- 
ческого самовозбуждения в системе (17) в плоскости «и-—\», как это принято при анализе уравнения Матье-Хилла 


[43]. На нижней иллюстрации рис. 2, 6 приведена аналогичная диаграмма для системы (8). 
Здесь =0/0%, Оу = [(<+71с()/ т? при анализе уравнения (6) ии=©/(%, О, = (с/т)? при анализе уравне- 


ния (8). Все приведенные диаграммы параметрического самовозбуждения относятся к системе, параметры которой 
соответствуют точке «1» на рис. 2, а. 
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а) 6) 
Рис. 2. Область устойчивости системы при и; , =9=0,1,5=1,2 
и диаграмма параметрического самовозбуждения 
Если рассмотреть эти же диаграммы для точки «2» на рис. 2, а, диаграмма параметрического самовозбужде- 


ния для уравнения (6) остается практически без изменения, а для уравнения (8) лепестки неустойчивого поведения 
существенно уменьшаются, как показано на рис. 3, а. Приведем также примеры диаграмм параметрического самовоз- 
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буждения системы для случая, когда нет угловых коэффициентов равных нулю. В этом случае параметрическое само- 
возбуждение за счет периодических изменений матрицы упругости и скоростных коэффициентов становятся взаимо- 


связанными. На рис. 3, 6 приведены диаграммы для системы (4), имеющей угловые коэффициенты: х, =0,6, > =0,8. 


На рис. 3, в — для угловых коэффициентов Х1 = 0,8, Х› = 0,6. 
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Рис. 3. Примеры диаграмм параметрического самовозбуждения системы фрезерования 


На рис. 3, г при Х, =0,8, Х. =0,6 учтены связи, формируемые в зоне контакта задней грани инструмента с 


обрабатываемой деталью, то есть для полной системы дифференциальных уравнений (3). 
Анализ результатов. По мере увеличения частоты вращения инструмента и приближения её к одной из соб- 
ственных частот динамической системы фрезерования наблюдаются явления параметрического самовозбуждения си- 


стемы. Если рассматривать динамическую систему (4) как систему с постоянными параметрами (50 =5® =9=0), 
то по мере увеличения скорости резания существует ограниченный диапазон, в котором стационарная траектория 
упругих деформационных смещений теряет устойчивость. Это обусловлено, как показано ранее [41], преобразованием 
симметричной составляющей матрицы скоростных коэффициентов математической модели динамической системы 
фрезерования из положительно определенной в отрицательно определенную. Это преобразование зависит, прежде 
всего, от запаздывания сил по отношению к деформационным смещениям и от крутизны уменьшения сил при возрас- 
тании скорости. В реальных системах, начиная с некоторой скорости резания, система вновь теряет устойчивость. При 
этом симметричная часть матрицы скоростных коэффициентов, как правило, остается положительно определенной. 
Как показано выше, это обусловлено параметрическим самовозбуждением системы, то есть в этом случае механизмы 
потери устойчивости имеют другую природу. Аналогичный эффект наблюдается и при уменьшении жесткости ин- 
струмента. Поэтому при выборе технологических режимов, прежде всего скорости резания, существует оптимальный 
диапазон, в котором система фрезерования остается устойчивой. 

Исследования показывают, что на условие самовозбуждения оказывают влияние практически все параметры 
динамической связи, технологические режимы (прежде всего скорость резания, влияющая на частоту вращения ин- 
струмента) и геометрия инструмента. Параметрическое самовозбуждение зависит от близости параметров системы к 


границе области устойчивости подсистемы при 6% =5 =9=0. Приближение параметров к этой границе не только 
распгиряет лепестки неустойчивого поведения системы, но и приводит к образованию этих лепестков в области малых 


Вестник Донского государственного технического университета 2016, №2(85), 32-42 





скоростей резания. К расширению области неустойчивости приводит и увеличение связанности двух взаимодейству- 
ющих через процесс резания колебательных контуров подсистемы инструмента. Наибольшей стабильностью стацио- 
нарные траектории обладают в системах, уравнение которых близко к уравнению Матье-Хилла. В этом случае первый 
лепесток области неустойчивости близок к двойной собственной частоте колебательного контура. Все перечисленные 
выше факторы снижают частоту вращения инструмента, начиная с которой наблюдается параметрическое самовоз- 
буждение. 

Несмотря на то, что дополнительная связь, формируемая в контакте задней грани инструмента с обрабатыва- 
емой деталью, имеет смысл нелинейной диссипации, за счет параметрического взаимодействия она также вызывает 
расширение области неустойчивости и смещение ее границы в область малых частот вращения инструмента. Это свя- 
зано с тем, что вариации скорости подачи, обусловленные периодическими упругими деформационными смещениями 
инструмента, за счет этой связи лишь увеличивают вариации параметров динамической системы. 

Приведенный материал образует информационную базу для построения систем динамической диагностики 
состояния процесса фрезерования. В отличие от точения, для которого предложены алгоритмы диагностирования [44], 
при фрезеровании эти алгоритмы должны быть дополнены оценками, основанными на стробоскопическом отображе- 
нии Пуанкаре с периодом, определяемом периодом повторения установившихся стационарных траекторий. 

Выводы. Параметрическое самовозбуждение динамической системы фрезерования, обусловленное периоди- 
ческими изменениями практически всех параметров динамической системы, являются одним из важных факторов 
потери устойчивости стационарных траекторий. Эти траектории фактически определяют параметры качества форми- 
руемой резанием поверхности. Приведенный материал позволяет наметить пути выбора технологических параметров 
и геометрии инструмента из условий устойчивости траекторий формообразующих движений инструмента относи- 
тельно заготовки при фрезеровании концевыми фрезами, тем самым обеспечивая проектируемое качество обрабаты- 
ваемой детали. Особенно актуален выбор параметров при создании программ ЧПУ обработки деталей сложной гео- 
метрической формы. 
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